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During the last 10 years, INTEGRAL made a unique contribution to the study of accreting mil- 
lisecond X-ray pulsars (AMXPs), discovering three of the 14 sources now known of this class. 
Besides increasing the number of known AMXPs, INTEGRAL also carried out observations of 
these objects above 20 keV, substantially advancing our understanding of their behaviour We 
present here a review of all the AMXPs observed with INTEGRAL and discuss the physical inter- 
pretation of their behaviour in the X-ray domain. We focus in particular on the lightcurve profile 
during outburst, as well as the timing, spectral, and thermonuclear type-I X-ray bursts properties. 
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1. Introduction 

The INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) was launched on 2002 
October 17 to observe extreme astrophysical events in the gamma-ray sky. During the first ten 
years in orbit, INTEGRAL contributed to the discovery of three accreting millisecond X-ray pulsars 
(AMXPs): IGR J00291-H5934 in 2004 ||2|, g], ICR J175 11-3057 in 2009 [|] , and IGR J17498- 
2921 in 2011 iQ]. INTEGRAL also contributed to the understanding of two recently discovered 
AMXPs: HETE 1900.1-2455 in 2005 [^] and Swift J1749.4-2807 in 2010 [jT^]. XTE J1807- 
294 was the first AMXP serendipitously observed with INTEGRAL IBIS/ISGRI between February 
and Mar-ch 2003 The AMXP SAX J1808.4-3658 was also obser-ved serendipitously from 
the beginning of its outburst in 2008 until back to quiescence (i.e., between September 18 and 
October 19, 2008). These INTEGRAL data still remain unpublished. The first outburst of the 
AMXP XTE J1751-305 was observed in April 2002 and lasted 12 days INTEGRAL observed 
the second outburst from this source in March 2005 the third in April 2007 [|], and the fourth 



in October 2009 ^. The outbursts from XTE J 175 1-305 were all relatively short and similar to 
each other (see [ ^3| ] and reference therein). 

At present, 14 AMXPs are known. All these sources are X-ray transients, i.e. they spend 
most of the time in a quiescent phase (X-ray luminosities of 10^^ - 10^^ erg s~^) and sometimes 
undergo X-ray outbursts reaching luminosities of 10^^ - 10^^ erg s"^. During these events coherent 



pulsations are observed with frequencies in the range between 182 and 599 Hz (see e.g. [25, 19, 18] 
for a review). Since the launch in 2002, INTEGRAL provided a contribution to our understanding 
of the AMXPs by performing at least one observation of these objects per year. In this contribution 
we review all the INTEGRAL results on the AMXPs achieved after ten years in orbit. 

2. The INTEGRAL imaging of AMXPs 

Thanks to the unprecedented INTEGRAL capabilities to image the hard X-ray sky at relatively 
high angular resolutions (12' for ISGRI and 3' for JEM-X), all AMXPs could be clearly identified 
in the mosaicked images at discovery and their X-ray emission unambiguously disentangled from 
that of other nearby hard X-ray emitters ^ 0, ^ 10 Hi- An example of mosaicked sky image 



around an AMXP in outburst is shown in Fig. 1 (left panel). 

These capabilities were particularly relevant for the observation of the source IGR J00291 -1-5934, 
which is located ~ 18' from the close-by persistent intermediate polar V609 Cas [^]. The AMXP 
source XTE J175 1-305 underwent a short period of enhanced X-ray activity during the decay of 
the outburst from the AMXP IGR J175 11-3057 located ~ 20' away from the former object In 
the case of the AMXP IGR J 17498-2921 INTEGRAL was able to disentangle the contributions of 
the sun^ounding objects and to clearly distinguish between the type-I X-ray bursts emitted by the 
source and those going off in the near-by bursters SLX 1744-300/299 and lA 1742-294 [|]. These 
sources are located at 0.9° and 0.86° from the IGR J17498-2921 position (see Fig. 1; left panel). 

3. The AMXP outburst lightcurves 

Most of the AMXPs displayed in the past very similar behaviours while in outburst. A repre- 
sentative Ughtcurve of an outburst from an AMXP is shown in Fig. 1 (right panel). The lightcurves 
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Figure 1: Left: INTEGRAL/lSGRl sky image of the field of view around IGR J17498-2921 in the 20-100 keV band 
for an effective exposure of 210 ks (see [|]). Right: RXTE/PCA (2-20 keV) outburst light curve of IGR J17511- 
3057. The diamonds correspond to observations in which both XTE J175 1-305 and IGR J1751 1-3057 were active. The 
arrows indicate the times of the detected X-ray bursts. The vertical dashed lines indicate the interval of the INTEGRAL 
observations (see ||9||). 



of these events are usually characterized by a fast rise (a few days at the most) and an exponential 
decay. The latter terminates with a break, after which a linear decay of the source X-ray flux is 
observed down to the quiescence level (see e.g., [jH, |6l, ^ |8], [T^, ^). These outburst profiles are 
commonly interpreted in terms of the disk instability model, taking into account the irradiation of 



the disk by the central X-ray source [|14|. This model has been applied to a sample of diff'erent 
LMXBs, including six AMXPs |22, ^,[1^, ^. In these works it was shown that the timescale of the 
decay light curve and its luminosity at a characteristic time are linked to the outer radius of the ac- 
cretion disk. The break ("knee") observed during the X-ray flux decay at the end of the outburst is 
believed to be a consequence of mass transfer onto the outer edge of the disk, as this supply of ma- 
terial to the compact object is effectively cut-off when the outer disk enters the cool low-viscosity 
state. The knee occurs at the lowest value of the X-ray luminosity at which the outer disk edge can 
still be kept hot by central illuminating source. The measured value of this luminosity can be used 
to estimate for the outer disk radius. A second independent estimate of this radius is given by the 
timescale of the exponential decay of the lightcui've. 



4. The AMXP X-ray spectra 



The available INTEGRAL JEM-X-fISGRI broad-band X-ray spectra of all AMXPs (XTE 
J1807-294, IGR J00291-F5934, HETE 1900.1-2455, IGR J17511-3057, Swift J1749.4-2807, and 



IGR J17498-2921), can be usually well fit by using a thermal Comptonization model compps [ |21[ ] 
in the slab geometry ||, ^, |8), 10, ||]. The main model parameters ai^e the Thomson optical 
depth tt ~ 1 - 2 across the slab, the electron temperature kT^ ~ 25 - 50 keV, the temperature 
kT\,\, ~ 0.5 - 1.0 keV of the soft-seed blackbody photons assumed to be injected from the bot- 
tom of the slab, the emission area Abb ~ 20 km^. In Fig. 2 (left panel) we show a representative 
AMXP unfolded spectrum and the corresponding best-fit compps model. 
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Figure 2: Left: Simultaneous INTEGRAL, XMM-Newton, and RXTE spectra of XTE J1807-294 fitted with an absorbed 
disk black body, diskbb, plus compps model (between March 20-22, 2003). The EPIC-pn and M0S2 spectra in the 0.5- 
10 keV range and PCA/HEXTE in the 3-200 keV and IBIS/ISGRI in the 20-200 keV range are shown. The diskbb 
model is shown by a dotted curve, the dashed curve represents the compps model, while the total spectrum is shown by a 
solid line. The lower panel presents the residuals from the best fit. See |^ for more details. Right: Pulse-profile collage 
of IGR J00291+5934 using data from RXTE/PCA (2-58 keV; panels a-e), RXTE/REXJE (20-151 keV; panels g-j), 
INTEGRAL/JEM-X (5-10 keV; panel f) and INTEGRALflSGRl (20-150 keV; panels k-n). Two cycles are shown for 
clarity. All profiles have high off-set values due to the intrinsic properties of the instruments. The error bars represent 1 
sigma statistical errors. All profiles reach their maximum near phase ~ 0.95. It is noteworthy the highly sinusoidal shape 
of the profiles for energies up to 100 keV. 

The spectra of AMXPs are very similar to those of the so-called "atoll" sources at low lu- 
minosities, where the X-rays are probably produced in the boundary/spreading layer near the NS 
equator (see e.g., [pT]]). These similarities can be explained by assuming that in both classes of 
sources the energy dissipation occurs in an optically thin medium (i.e. the accretion shock and 
boundary/spreading layer) and the spectral properties are determined solely by energy balance and 
feedback from the NS surface (which in turns provides cooling in the form of soft photons). 

5. The AMXP timing above 20 keV 

In recent years, INTEGRAL proved to be crucial to perform timing analysis of the high-energy 
emission from AMXPs in the hard X-ray domain (>20 keV). In the case of IGR 100291-^5934, 
IGR J17511-3057, and IGR J17498-2921, INTEGRALflSGRl was capable of measuring the pulse 
profile at frequencies of 1.67 ms, 2.5 ms, and 4.08 ms, respectively up to energies of 150keV 
[§, ^. In the case of IGR J00291-I-5934 see Fig. 2 (right panel). These measurements per- 
mitted to perform an analysis of the pulsed fraction up to high X-ray energies and to confirm for 
the first time the increase of the pulsed fraction with energy in an accretion powered AMXP [^]. 
The latter could be succesfully interpreted in teniis of Doppler effects on the exponentially cut- 
off Comptonized photons produced from the hot spot on the neutron star surface (see Fig. 3; left 
panel). Soft phase/time lags (low-energy pulses lag behind the high-energy pulses) with a complex 
energy dependence could also be measured and linked to the same effect. According to previous 
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Figure 3: Left panel: The pulsed fraction of IGR J00291+5934 using pulsed/total flux measurements from 
PCA/HEXTE and JEM-X/ISGRI (open circles, PCA pulsed flux relative to JEM-X total flux; open squares, HEXTE 
pulsed flux relative to ISGRI total flux; filled squares, ISGRI pulsed flux relative to ISGRI total flux). The pulsed frac- 
tion gradually increases from ~ 6% at ~ 6 keV to ~ 12-20% at - 100 keV [|]. Right panel: IGR J0029 1+5934 time lags 
as a function of energy in the 2 - 100 keV energy range combining RXTEfPCA (2-35 keV, open circles), RXTEfHEXTE 
(20-60 keV, filled diamonds), and INTEGRAL/lSGRl (20 - 100 keV, filled squares) measurements.The time lags first 
increase and then decrease slightly, saturating above 15 keV, and possibly reaching zero around 50 keV [p|]. 



measurements performed on the X-ray emission from SAX J1808.4-3658 and XTE J1751-350, 
the time lags revealed by INTEGRAL monotonically increased with energy and saturated at about 
10-20 keV. The lags are most probably due to the different emission patterns of the blackbody and 
Comptonization components [20, n], |l3]] combined with the action of the Doppler effect. This in- 
terpretation is supported by the energy dependence of the lags, which grows until the contribution 
of the blackbody becomes negligible. At higher energies, the soft lags can evolve further. This is 
due to the fact that higher-energy photons suffer more scatterings and the latter produce notable 
variations in the emission pattern as a function of the energy (see e.g, ^ and reference therein). 
An example of measured time lags in the AMXP IGR J00291 -1-5934 is shown in Figure 3 (right 
panel). In this case, the time lags first increase and then decrease slightly, saturating above 15 keV, 
and possibly reaching zero around 50 keV (see also 

By using RXTE/PCA data for IGR J0029 1-1-5934 in outburst, it was possible to measure for 
the first time a decrease of the NS spin period while the source was accreting material [^. This 
provided an evidence in favor of the so-called "recycling" scenario, according to which old neutron 
stars in LMXB become millisecond X-ray pulsars when they are spun-up to such periods by the 
transfer of angular momentum from the accreting matter. This result was broadly communicated to 
scientific community through an ESA press release and global coverage in the popular press (see 
http://www.esa.int/esaSC/SEMWSAA5QCE_index_0.html). 

The complementary INTEGRAL and RXTE observations provide an opportunity to cross- 
calibrate the absolute timing of the two instruments. The cross-correlation of the folded solai^ 
system baiycentered event lists collected during the outburst of IGR J175 1 1-3057 from ISGRI (15- 
90 keV, 60 bins), RXTE/FCA (15-60 keV, 60 bins) and RXTE/HEXTE (15-90 keV, 60 bins) in the 
same energy band gave time shifts of -49 +/- 114 /us, and -69 +/- 49 /is, respectively [|T5[ ||, |, ^. 
INTEGRAL event times seemed to be 50 fis ahead with respect to the RXTE ones. Note, that 
this is almost the value of 47 yus that has been added in the process of updating/regenerating the 
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INTEGRAL time correlation files and used also to update the ISGRI User Manual. So, from an 
experimental point of view there is no justification to apply this shift of 47 /us in the INTEGRAL 
con^elation files, but within eiTors everything still matches. We can state that INTEGRAL ISGRI and 
RXTE PCA/HEXTE timing measurements ai^e consistent in absolute sense at the 50 fis level. We 
have thus the absolute timing information to perform a confident timing analysis. However, only 
with the RXTE monitoring observations over long time spans (to establish orbital parameters) an 
ephemeris could be constructed and used for the INTEGRAL/lSGRl events data. A blind searches 
with INTEGRAL will not lead to establish the AMXP nature. One long XMM-Newton observation 
may establish the AMXP nature, but can not yield an ephemeris without a long monitoring pro- 
gram. Nowadays without monitoring observations and without sensitive X-ray instruments with 
ms-timing capabilities the discovery of new AMXPs is unlikely. 



6. The AMXP observed X-ray bursts 
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Figure 4: Left panel: The JEM-X (3-20 keV, upper panel) and IBIS/ISGRI (18^0 keV; lower panel) background 
subtracted lightcurve of a bright X-ray burst detected from HETE J1900. 1-2455. The time bin is 0.5 s. At high-energy 
the burst shows strong evidence of photospheric radius expansion [Q]. Right panel: Triangles in the figure represent the 
observed burst recurrence times shown as a function of local mass accretion rate (the corresponding flux is reported in 
the upper x-axis) for IGR J1751 1-3057. We also show the best-fit power-law model. The recurrence time increased with 
time roughly as (Fpe,,s,bol)"' ' 

INTEGRAL observed several type-I X-ray bursts from the sources HETE 1900.1-2455, Swift 
J1749.4_2807, IGR J175 11-3057, and IGR J17498-2921 [0, [l^, |, |]. Themionuclear (type-I) 
X-ray bursts ai^e produced by unstable burning of accreted matter on the NS surface. The spectrum 
from a few keV to higher energies can usually be well described by using a blackbody model with 
temperature kTb\, wl-3 keV. The energy-dependent decay time of these bursts is attributed to the 
cooling of the NS photosphere and results in a gradual softening of the burst spectrum with time 
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(see e.g., [ |T^ for a review). Long (~220 ks) uninterrupted INTEGRAL observations of AMXPs 
in outbursts permitted to detect series of type-I bursts and constrain accurately their recunence 
time. In addition, the dependency of the recunent time on accretion rate and the ignition depth 
could be measured by making use of combined observations caiiied out by INTEGRAL, RXTE, 
Swift, Chandra, and XMM-Newton. In all the above mentioned sources, the short burst profiles 
provided indication of hydrogen-poor material at ignition, suggesting either that the accreted ma- 
terial is hydrogen-deficient, or that the CNO metallicity is slightly higher than the solar value in 
the AMXPs. For instance, for IGR J1751 1-3057 (see Fig. 4 right panel) the variation of the burst 
recurrence time as a function of m (inferred from the X-ray flux) is much smaller than predicted by 
helium-ignition models [^]. 

7. Summary 

INTEGRAL was able to observe at least one AMXP per year and discovered three objects be- 
longing to this class. By combining high-energy INTEGRAL data with those obtained from RXTE, 
XMM-Newton, and Swift, it was possible to study the timing and bursts properties of the AMXPs 
in a wide energy range (0.3-150 keV). Our reseai^ch group has been leading INTEGRAL observa- 
tional campaign of AMXPs since the early operational phases. Our AMXP Target of Opportunity 
(ToO) proposals were approved during AO2-A10 cycles, and during the past few years our ToOs 
were triggered five times. We obsei-ved IGR J00291-I-5934 in 2004, HETE 1900.1-2455 in 2005, 
IGR J175 11-3057 in 2009, IGR J17498-2921 in 2011, and Swift J1749.4-2807 in 2010. These 
observations were possible thanks to the prompt reactions to our requests from C. Winkler and 
the INTEGRAL staff, to whom we are grateful. We are also grateful to all TAG members who 
positively evaluated our research programs over years and permitted crucial advancements in the 
understanding of the physics of AMXPs. 
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